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RESEARCH SUMMARY 

This paper documents a computer model designed to 
simulate stand development in stands affected by the 
Douglas-fir tussock moth. The simulation model is actu- 
ally a combination of two independently developed 
models: the Stand Prognosis Model and the Dougjas-fir 
Tussock Moth Outbreak Model. This Combined Model 
can be used to assess the likely consequences of both silvi- 
cultural treatments and tussock moth control activities for 
stands in the Northern Rocky Mountains, using existing 
forest inventories. It can be used as a tool for long-range 
timber management planning because it displays the pro- 
jected results of alternative strategies for the management 
of forests affected by the tussock moth. This integrated 
approach permits direct comparisons of various manage- 
ment and tussock moth control strategies in terms of stand 
volume development over time, rather than an inter- 
mediate effect such as defoliation. The flexibility of the 
model has also proved valuable in examining the impor- 
tance and sensitivity of various assumptions in the Com- 
bined Model, and thus is useful in pointing out future re- 
search needs. 

This paper covers four major areas: (1) an overview and 
brief discussion of the Combined Stand Prognosis and 
DFTM Outbreak Model is given; (2) a description of the 
information needed to use the Combined Model is given, 
which includes documentation and discussion of the input 
options; (3) the output and information produced by the 
Combined Model is discussed; and (4) numerous examples 
are presented that illustrate the behavior and sensitivity 
of the Combined Model when major input options are 
varied. 
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INTRODUCTION 


Purpose 


This paper reports on a computer model designed to simulate stand development in 
forest stands affected by the Douglas-fir tussock moth (DFTM), Orgyia pseudotsugata 
(McDunnough), which is a defoliator of true firs, Abies spp., and inland Douglas-fir, 
Pseudotsuga menziesii var. glauca (Beissn.) Franco. The simulation model is actually a 
combination of two independently developed models: the Stand Prognosis Model (Stage 
1973; Wykoff et al. 1982), and the DFTM Outbreak Model (Overton and Colbert 1976 et 
seq.; Colbert et al. 1979, 198T; Overton et al. 198U). The Stand Prognosis Model was 
developed by the Intermountain Forest and Range Experiment Station in Moscow, Idaho. 
The DFTM Outbreak Model was developed jointly by Oregon State University and the 
Pacific Northwest Forest and Range Experiment Station in Corvallis, Oreg. (see Colbert 
[1978] for a short history of the development of the Outbreak Model). The effort that 
resulted in the combining of these models was sponsored by the Expanded Douglas-fir 
Tussock Moth Research and Development Program. 


tion of the program options and a description of propam input, 

• to discuss the output and information produced by the Combined Model, and 

• to provide examples that illustrate the behavior and sensitivity of the Combined Model 

when major input options are varied. 


•r'f^lhpr^ 1 J W S Overton and C. White. 1981 . Behavior of the Douglas-fir tussock moth outbreak popu- 
tion^ 

W S B E Wickman, and R. R. Mason. 1981 . Nature, organization, and content of “ 
jpulallon ieaks of'the Douglas-f.r lussock moth. IM p, Manuscript in process and on nie at Pac. Northweri 
ar. and Range Exp . Stn . , Corvallis , Oregon . 
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OVERVIEW OF THE COMBINED MODEL 
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The Stand Prognosis 
Model 


The Douglas-fir 
Tussock Moth (DFTM) 
Outbreak Model 


Another advantage of such a simulation model is that the user can easily see the effect of 
varying most assumptions that have been made in describing stand and insect conditions. 
Those assumptions that are supported by the least reliable information can be examined in 
greater detail; an ad hoc sensitivity analysis can usually be performed rather easily and 
should point out the assumptions that are critical, as well as those that are not. If necessary, 
field sampling can then be used to increase the reliability of those assumptions that appear 
to be either critical or weak. 


Stage (1973) described a Stand Prognosis Model that has become an increasingly useful 
tool in examining long-term stand management alternatives (Stage et al. 1980). The Stand 
Prognosis program is an individual tree-based stand model designed to simulate the de- 
velopment of the mixed-species even- and uneven-aged stands commonly found in the 
Northern Rocky Mountains. The simulator consists of separate component models for tree 
diameter growth, height growth, crown ratio, and mortality. The projection is produced by 
the repeated addition of periodic increments on diameter, height, and crown ratio to the 
initial dimensions of the inventoried trees. Numbers of trees are reduced by periodic mortal- 
ity rates. Normal period length is 10 years. 

Usual input for the Prognosis program is a list of tree and stand characteristics obtained 
from a standard stand inventory (Stage and Alley 1972). The diameter and species of each 
sample tree is required, while only a subsample of total height, crown ratio, and past gro^h 
is anticipated (but not required); spatial information describing individual tree locations is 
not used. Site characteristics sampled include slope, aspect, elevation, habitat type, and 
geographic location; measures of stand density are calculated internally from the list of tree 
characteristics. Information describing the sampling design is also required so that the 
number of trees per unit area represented by each sample tree can be calculated (Stage 
1978d). Stand statistics, such as basal area or volume per unit area, are then calculated by 
summing the corresponding tree characteristics, weighted by the number of trees repre- 
sented by each tree. Thus, the Prognosis Model can be used to display relevant stand 
statistics versus time even though the basic unit in. the model is an individual tree. A more 
detailed discussion of the use of the Stand Prognosis Model in timber management applica- 
tions is given by Stage (1978 et seq.) and Stage et al. (1980). 


The DFTM Outbreak Model (Overton and Colbert 1976 et seq.; Colbert et al. 1979, 198) 
[see footnote U; Overton et al. 1981 [see footnote 2]) simulates the course of events during 
a tussock moth outbreak on a collection of 1000 in* midcrown sample branches (see Mason 
1970) The outbreak is assumed to be 4 years in duration, with each year corresponding to a 
distinct phase in the outbreak (Mason and Luck 1978). The model was calibrated using datj 
obtained during the last (1971-74) Blue Mountains outbreak in northeastern Oregon (see ^ 
Mason 1976, 1978; Wickman 1978 et seq.; Beckwith 1978). The major processes considered 
are insect survival, growth, and feeding and the associated host defoliation; an 
redistribution of insects between sample branches is also considered. To run the D^M 
Outbreak Model, the following information is needed: insect population denpty, biomass, 
of the new foliage and either percentage new foliage or biomass of the old foliage; the ^ 
foliage information is needed for each of the host species considered (Douglas-fir and grap< 
fir'i Colbert and Wong (1979) have detailed the procedures for running this simulation 
iS^ly sland Prognosis Model; see » et al. (19^ 
documentation. Additional discussion has been provided by Overton and Colbert (19J8a.i: 
b c) and Colbert (1978). The version of the DFTM Outbreak Model (version 3.1) us^^| 
subroutines in the Combined Model discussed in this paper is the same as that descnbed by 
Colbert and Wong (1979) and Colbert et al. (1979). except that the output tables have beejt 

deleted. I 
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The Combined Model 


Tree Defoliacion Effects 


To properly use these two rather disparate models in conjunction, an understanding of 
the following topics is necessary: 

1 . Tree defoliation effects considered by the model, 

2. Foliage biomass classification model, 

3. Tree class compression, 

4. Methods for allocating first instar larvae to tree classes, 

5. Probability of outbreak model, 

6. Outbreak control and salvage options available. 

Each of these topics was discussed in detail by Monserud {1978a). A brief discussion will 
be given here. 


A tussock moth outbreak can affect normal tree development in three major ways: 

1. Growth in height and diameter can be retarded, 

2. Total height and volume can be reduced because of top-kill, and 

3. Probability of mortality can increase. 

Although there may be other effects (for example, fertilization due to rapid nutrient turn- 
over following defoliation), only these three are simulated in this model. The research basis 
for the quantification of these defoliation effects is due almost entirely to the work of B. E. 
Wickman, at the Corvallis Forestry Sciences Laboratory (for examples, see Wick man 
1978a,b; Wickman et al. 1980). 

Output from the DFTM Outbreak Model consists of percentage midcrown branch de- 
foliation in Phase II and Phase III of the outbreak for each tree or class of similar trees. 
Maximum defoliation on the midcrown branch is then converted to percentage tree defolia- 
tion by the function described by Overton and Colbert (1978b), which is graphed in figure 2. 
Two important characteristics of this function are: 

1 . Percentage of tree defoliation remains essentially zero until midcrown branch defolia- 
tion exceeds 55 percent; 

2. Percentage of tree defoliation then increases rapidly until complete defoliation of 
both branch and tree is reached. 



Figure 2 -Percentap tree crown defoliation versus percentage maximum midcrown 
sample branch defoliation (from Overton and Coibert 1978b). 



Foliage Biomass 


Tree Class Compression 


Allocation of First Instar 
Larvae 


it is very important to be aware of the behavior of this function when experimenting with 
the Combined Model, for a seemingly minor change in initial conditions could result in a 
dramatic change in tree defoliation if branch defoliation is between 70 and 90 percent. And 
any combination of control options or initial conditions that keep branch defoliation from 
exceeding 55 percent will produce results essentially identical to a simulation with no tus- 
sock moth outbreak at all. 

Percentage of tree defoliation for a particular tree or class of similar trees is then used to 
index the table of defoliation effects (see Table 1 1-2 on p . 226 of Colbert and Campbell 
1979, and discussion by Monserud 1978a, p. 65-66); note that the sampling basis for the 
defoliation effects is approximately 5 years (Wickman 1978a). This table contains the fol- 
lowing species-specific information for seven tree defoliation classes: 

1 . Probability of direct mortality caused by tussock moth; 

2. Probability of secondary and background mortality, 

3. Probability of top-kill occurring on surviving trees for each of five top-kill classes; 

4. Percentage of diameter and height growth reduction. 

With the exception of calculating volume loss due to top-kill (Monserud 1980, 1981) the 
modification of tree characteristics resulting from defoliation is straightforward. 


The DFTM Outbreak Model of Overton and Colbert (1976) simulates the course of 
events during an outbreak on a collection of midcrown sample branches of 10(X)-inS as 
described by Mason (1970). The sample branch has only three essential characteristics: hosi 
species (Douglas-fir or grand fir), percentage new foliage, and total foliage biomass. Since 
the basic unit in the Stand Prognosis Model is a tree, and the basic unit in the Outbreak 
Mode! is a branch, a linkage between these units is needed to combine the two models. 
Species-specific equations for predicting foliage biomass characteristics of the .sample 
branch, from tree and stand information available to the stand model, facilitate this link- 
age. Available equations were developed by C. R. Hatch at the University of Idaho (Hatch 
and Mika 1978; Monserud 1978a, p. 49-56). Unfortunately, the data base supporting these 
equations is quite limited, and should be improved. By sampling for the mean and standarc 
deviation of the necessary foliage characteristics, the user is likely to increase the accuracy 
of the resulting simulations. 


The current version of the Stand Prognosis Model (Inland Empire version 4.0, as docu- 
mented by Wykoff et al. 1982) can handle up to 1 ,350 individual tree records. Even when 
the number of records in a sample stand is small, the record tripling logic described by Stag 
(1973, p. 12-13) Will usually result in several hundred tree records being projected. The 
DFTM Outbreak Model, however, only views trees in two dimensions (foliage biomass and 
percent new foliage) for each host species. Considerable computer time can therefore be 
saved by “compressing” trees with similar foliage biomass characteristics into the same tree 
cla^s before the Outbreak Model is called (see discussion preceding the NUMCLASS key- 
word in a following section). Each “tree class” represents one or more tree records in the 
stand model and is represented by one midcrown sample branch (Mason 1970) in the Out- 
break Model. The current version of the Combined Stand Prognosis and DFTM Outbreak 
Model allows for a maximum of l(X)tree classes; this arbitrary limit appears to provide quit 
adequate resolution for simulating defoliation effects. ! 


In addition to the foliage biomass information, the Tussock Moth Outbreak Model re- i 
quires the specification of the number of established first instar larvae for each tree class I 
at the start of an outbreak (Phase I). A tree class is simply a class or group of similar trees i 
represented by a single midcrown sample branch . Note that ‘ ‘ first instar larvae* ‘ in this | 
paper should be considered synonymous with “viable eggs’ * in earlier papers (e.g., Monsej 
rud 1978a, b). i 


Probability of Outbreak 


Two methods—assumptions, actually— are currently available for allocating levels of first 
instar larvae to the tree classes: the first is random and the second is deterministic. These 
methods allow for examining the various assumptions concerning the between-tree distribu- 
tion of larvae in a stand. 

With the random method, the number of larvae assigned to a particular tree class is 
drawn from a random normal distribution with specified mean and standard deviation. This 
option is quite useful when a sample of early instar larvae is available for a stand, assuming 
that no relationship can be found between larval density and observable tree characteristics. 
The random method also contains a feature that allows the user to vary the mean larval den- 
sity if multiple outbreaks are being simulated. 

The deterministic allocation method assigns three specified larval levels or densities to 
each ordered third of the tree classes, after the tree classes have been sorted by average 
diameter; the first tree class to have larvae assigned to it has the largest average diameter. 

Wickman (1978a) reported on the distribution of mortality by diameter class in the last 
Blue Mountains outbreak. Either of these larval allocation assumptions could be used to 
mimic the evenly distributed pattern (with respect to diameter) of mortality for grand fir. 

The deterministic assumption, however, would be most amenable to reproducing the distri- 
bution of mortality reported for Douglas-fir, which was more concentrated in smaller diam- 
eter trees. 

The actual levels of larvae specified with any of these options should be based on the 
most reliable estimates available, for this is the most important variable in the outbreak 
model. The choice of the assumption or method most appropriate for a particular situation 
should— if possible— be determined by analyzing tussock moth inventory data in relation 
to tree diameter or size. If this is not possible, we recommend using the random allocation 
method. 


The DFTM Outbreak Model is just that — it simulates the course of events during a 
tussock moth outbreak. It was not designed to model population dynamics when the 
tussock moth is not at outbreak levels. Thus, the crucial decision of whether or not to in- 
voke an outbreak must be made before the DFTM Outbreak Model is called. This forces 
consideration of the probability of an outbreak occurring in a given stand, in a given year. 

The salient features of past outbreak patterns (Stage 1978a) are; outbreaks appear to 
be synchronized over large areas; intervals between outbreaks are usually at least 8 years; 
some stands are involved in repeated outbreaks, while others are involved only once; 
and not all the stands with similar conditions and histories are involved in a given out- 
break. 

These four features can be represented by the two-step process described by Stage 
(1978a). First, a sequence of dates is stochastically determined that represent the times when 
large-area outbreaks are to be simulated. This sequence represents the temporal probability 
of outbreak. Second, the relative probability with which a stand of particular attributes can 
be expected to show defoliation at the time specified by step one is then determined; this is 
the spatial probability of outbreak, conditional on the temporal probability of outbreak. 
The model used to predict the spatial probability of outbreak was developed by Heller et al. 
(1977), and uses both physiographic variables (slope, aspect, elevation, topographic posi- 
tion) and stand variables (crown closure, percent host species, average crown diameter). A 
similar outbreak mode! has recently been developed for the Palouse Ranger District (Clear- 
water National Forest) northeast of Moscow by P. B. Mika and J. Moore (personal com- 
munication 1979; see Stoszek et al. 1981 for a related analysis) at the University of Idaho, 
Mika and Moore’s model was developed from data collected on ground plots, whereas data 
used by Heller et al. (1977) was obtained from aerial photographs in the Blue Mountains 
and Col^^e areas. The long-term accuracy of these models for predicting the conditional 
probability that a given stand will be involved in a regional outbreak is unfortunately 
unknown, for they are developed from data collected in only one regional outbreak 
(1971-74). Such models are also likely to be conditional on the specific (and unquantified) 

climatic factors associated with the 1971-74 outbreak. 



Outbreak Control and 
Stand Management Options 


The use of the probability-of-outbreak model is not a necessary feature of the combined^ 
system; an option is available for simulating any predetermined sequence of outbreaks. This 
option is quite useful for making retrospective comparisons of management alternatives 
when a stand's outbreak history is known. 

The DFTM Outbreak Model can be used to simulate a number of control options— 
biological as well as chemical— by altering mortality rates at specific occasions in the out- 
break. In the combined model, a simulated chemical control can be applied to any instar in 
any phase of the outbreak with any efficacy. Biological control options are available for ap- 
plying nuclear polyhedrosis virus (NPV) in either Phase II or Phase III of the outbreak. 

A wide variety of silvicultural options are also available in the Stand Prognosis Model. 
Thinning is simulated by reducing the number of trees per acre represented by the tree 
records until a user-specified thinning target is reached. The thinning target can be specified 
as: 

a residual number of trees per acre; 
a residual basal area per acre; 
a segment of the diameter distribution; 
a percentage of full stocking; 

a prescription where specific tree records are coded for cutting. 

The first three targets can be reached by thinning either from above or below. These thin- 
ning options can be implemented in any cycle of a simulation. Options are also available for 
specifying species preferences for harvesting that would allow for selectively removing host 
species. Use of the thinning options is discussed by Wykoff et al. (1982), 

An additional harvesting option is also available in the Combined Model: a salvage thin- 
ning operation immediately following the outbreak can be requested . All host trees that _ 
have been defoliated more than a (user-supplied) minimum percentage tree defoliation will 
be salvaged. Thus some of the volume that normally would be lost due to tussock moth can 
be recovered. Of course, this will usually result in harvesting addition^ trees that are par- 
tially defoliated, but not killed— as happened in the last Blue Mountains outbreak. 


DOCUMENTATION OF INPUT OPTIONS 


The Combined Stand Prognosis/DFTM Outbreak Model utilizes a keyword system for 
specifying program options. The keywords are intended to simplify the process of specify- 
ing the various features of the model that will be implemented in a given simulation or run. 
The keyword system is analogous to a high level (albeit simple) computer language in which 
each instruction or keyword is translated into a number of complicated instructions. To in- 
voke a particular option, the user simply inserts a short keyword in the runstream that is 
sent to the computer. For example, the keyword NPV2 specifies the application of nuclear 
polyhedrosis virus in Phase II (2) of an outbreak; although four separate tussock moth mor- 
tality rates need to be altered to simulate this option, these rates do not need to be provided 
by the user. The user is required to supply much less information to effect the desired result. 
Furthermore, the order that the information or keywords are submitted is usually not im- 
portant; of course there are exceptions to this rule (for example, the DFTM keyword must 
precede all other tussock moth keywords). This keyword system is used extensively in the 
Combined Model, and has greatly simplified the process of preparing a simulation. 

An additional simplifying feature of the keyword system is that default values exist for 
almost all keywords. Only keywords for nonstandard options need be specified, and the 
only numeric (parameter) values necessary on such keyword cards are those that differ from 
the default parameter values. In addition, all options and parameter values are reset to the 
default values after projecting a given stand in a multi-stand simulation. 



Rules for coding keywords: 

1 . All keywords start in column 1 of the keyword card or record. 

2. Numerical values (termed “parameters”) needed to implement an option are con- 
tained in seven numeric “fields” that are each 10 columns wide, beginning in column 
1 1. A decimal point should be punched for all numeric values that are not integers, 
and integer values should either be right-justified in the numeric field or followed by a 
decimal point. If a decimal point is provided, the actual location of the numeric value 
within the 10 column wide field is unimportant. 

3. Blanks that are coded in the numeric fields are not treated as zeroes. If a blank field 
is found, the default value will be used. If zeroes are to be specified, they must be 
punched. Thus, only the numeric values that are different from the default parameter 
values need be specified (in the appropriate field, of course). 

4. When two or more conflicting options are encountered, the last one specified will be 
used. 

5. The first tussock moth keyword must be DFTM, and the last tussock moth keyword 
must be END. 

This paper documents only those keywords used to implement the various features of the 
tussock moth portion of the Combined Stand Prognosis/DFTM Outbreak Model. Docu- 
mentation for the keywords used by the Stand Prognosis Model is provided by Wykoff et 
al. (1982). Although these additional keywords are necessary to run the Combined Model, 
they are not listed here to avoid duplication with Wykoff et al. (1982). 

Keywords will be grouped into the following five categories in the subsequent sections: 
program execution options, outbreak timing options, outbreak initial conditions options, 
tree cl^s compression and redistribution options, and outbreak control options. Each sec- 
-tion will contain a definition and short description of relevant keywords. Most sections will 
end with examples illustrating the use of keywords to simulate specific situations. 


Program Execution 
Options 


The keywords for controlling program execution options serve four general functions. 
The DFTM and END keywords signal the Stand Prognosis Model that tussock moth 
keywords will follow and have ended, respectively. NODFRUN and NOG F RUN can 
be used to exclude either Douglas-fir or grand fir as a host for tussock moth. DEBUG, 
DEBUTREE, PUNCH, REPORT, and DATELIST provide supplemental (and occasionally 
voluminous) output useful for examining the program’s behavior in greater detail than that 
afforded by the usual output. And RANNSEED allows the user to choose a different ran- 
dom number sequence in any of the routines related to tussock moth. With this last option, 
a user can assess the magnitude of the variability associated with the stochastic components 
of the Combined Model. We anticipate that the only keywords in this section that most 
users will need to use are DFTM, END, and RANNSEED. 


Keyword Keyword Description 


DFTM 

END 


NODFRUN 

and 

NOGFRUN 


Signal the Stand Prognosis Model that tussock moth keywords follow; 
use of this keyword is mandatory. 

Signal the Combined Model that the preceding group of tussock moth 
keywords has ended. Any number of groups of tussock moth keywords 
can be used, provided each group begins with the DFTM keyword and 
ends with the END keyword. Use of this keyword is mandatory. 

These keywords inhibit the DFTM Outbreak Model from simulating 
the activity of tussock moth on Douglas-fir (NODFRUN) and/or grand 
fir(NOGFRUN). 



DEBUG 

DEBUTREE 

PUNCH 


Field 1: 


REPORT 


Field 1: 


DATELIST 

RANNSEED 


Field 1: 
Field 2: 
Field 3: 


A large amount of intermediate output detailing the operation of 
the Combined Modei will be printed. This and the following option 
should rarely be needed by the normal user. 

Long tables of intermediate values associated with individual tree 
records will be printed. 

The input values and parameter arrays needed by the DFTM Outbreak 
Model as well as the defoliation levels by tree class, species average, and 
stand average are written in card image format to a separate output 
unit. 

The FORTRAN data set reference number where the supplemental 
output is written; there is no default.The user must specify a valid 
number and the corresponding job control statement. 

Control the amount of tussock moth output generated by the Com- 
bined Model; default is 2. 

The report level, where: 

0 = no tussock moth output will be generated; 

1 = only the DFTM Outbreak Summary Table will be printed; 

2 = All normal tussock moth output tables will be printed (de- 
scribed in the INFORMATION PRODUCED section). 

All tussock moth subprograms (i.e., subroutines and functions) 
and common blocks are listed with the date each was most recently 
revised. 

The pseudorandom number generator (Marsaglia and Bray 1968) used 
by the Combined Model has three seeds. These seeds initialize the ran- 
dom number generator, and are set at the beginning of each simulation. 
As a result, the random numbers will be generated in the same order 
each time a runstream is submitted. Consequently, identical projections 
of a single stand made in separate runstreams will have identical results. 
You can introduce some random variation by replacing one or more of 
the seeds. Since the new seeds should be odd integers, 1 will be added 
to any even numbers that are used to reseed the random number gener- 
ator. Note that this random number generator is used only by the tus- 
sock moth related subroutines in the Combined Model; the Stand 
Prognosis Model has a separate, but identical, random number gener- 
ator that is unaffected by this reseeding. 

First seed; the default is 1409859205. 

Second seed; the default is 402656419. 

Third seed; the default is - 328609067 . 


Outbreak Timing 
Options 


As discussed in the “Probability of Outbreak” section, specifying the time periods that 
the DFTM Outbreak Model is to be called is a two-step process. The first step schedules the 
occurrence of regional tussock moth outbreaks. This step can be handled in two different 
ways: manually (deterministicaUy), using the MANSCHED keyword, or randomly, using 
the RAN SC H ED keyword. The second step selects which regional outbreaks will include 
the subject stand. The M ANSTART keyword specifies that the subject stand will be includ- 
ed in aU regional outbreaks; the RANSTART keyword makes this determination random, 
conditional upon the stand’s probability of outbreak. The PROBMETH keyword specifies 
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the method for calculating the stand’s conditional probability of outbreak. The TOPO and 
ASH DEPTH keywords provide information on topographic position and ash depth, re- 
quired by some of the options available with the PROBMETH keyword. 


MANSCHED 

Field 1; 

RANSCHED 

Field 1; 

Field 2: 

Field 3: 
MANSTART 

RANSTART 

PROBMETH 

Field): 


Manually specify either the calendar year or the cycle number in which 
a regional outbreak will occur; the default is for NO regional outbreak.? 
to occur. If more than one regional outbreak is to be scheduled, use ad- 
ditional MANSCHED keywords. 

Either the year or the cycle number in which a regional outbreak will 
occur; default is for no regional outbreak to occur. NOTE: The 
Combined Model assumes that a number in Field 1 is a cycle number 
if it is less than or equal to 40— the maximum number of cycles (i.e., 
growth projection periods) allowed by the Stand Prognosis Model. 

Invoke the random automatic scheduling process which will sto- 
chastically generate a list of regional outbreaks which occur during the 
simulation period (see Stage 1978a). This is done by drawing from a 
random Bernoulli process with a specific minimum waiting time (Field 
1) and a specific event probability (Field 2); the process begins with the 
year of the last regional outbreak (Field 3). 

The minimum waiting time between regional outbreaks; default is 30 
years. 

The event probability used in the random Bernoulli process; default 
is 0. 1. This is essentially the annual probability of a regional out- 
break given that the minimum waiting time since the last outbreak 
has been exceeded. Note that the expected value of the time (T) bet- 
ween outbreaks is T = M + (1/P) - 1 , where M is the minimum 
waiting time (Field 1) and P is the event probability (Field 2). 

The calendar year of the last regional outbreak; default is year 1492. 

Specify that the DFTM Outbreak Model will be called whenever a 
regional outbreak is scheduled. MANSTART is the default “start” op- 
tion. 

Stochastically determine if the subject stand will be included in the 
regional outbreak. This is done by calculating a conditional probability 
(determined by the method specified in Field 1 ofthe PROBMETH 
keyword), that the subject stand will be infested by tussock moth, given 
that there is a regional outbreak. If this conditional probability is 
greater than a uniform random number between 0 and I , then the 
regional outbreak includes the subject stand. 

Select a method for calculating the conditional probability that the 
subject stand will be included in a regional outbreak. This keyword is 
normaUy used in conjunction with the RANSTART keyword. If this 
ke^ord is used with the MANSTART keyword, the conditional prob- 
ability of outbreak calculations will be made and printed in the “DFTM 
Outbreak Summary Table” (discussed in a later section), but wiU other- 
wise be ignored by the program. 

The conditional probability calculation method (default is method 
1), where: 

1 = Use the model developed by Heller et al. (1977), which is a 



function of elevation, slope, aspect, topographic position, stand 
closure, proportion of stand in host, and average crown width. The 
last of these three variables are calculated as functions of crown 
competition factor. This mode! was calibrated using aerial photo in- 
terpretation data from the Blue Mountains of northeastern Oregon. 
See the TOPO keyword below for details concerning the 
topographic position specification. 

2 = Use a model developed by Mika and Moore (personal com- 
munication, 1979) which is a function of topographic position (see 
TOPO), ash depth in inches (see ASH DEPTH), total basal area, and 
proportion of the stand in grand fir. This model was calibrated using 
data collected from the Palouse Ranger District of the Clearwater 
National Forest, northern Idaho. 

3 = Use a model similar to 2 above, except that ash depth is not 
used. This model was developed using the same data used to develop 
model 2. 

Field 2: The conditional probability scaling factor.^ The default scaling fac- 

tor is 1 .0; a value of 0.5 would reduce the calculated conditional 
probability by half, and a value of 2.0 would double the calculated 
value. 

TOPO This keyword is used in conjunction with the PROBMETH keyword 

to enter the numeric code specifying the topographic position of the 
stand. When the conditional probability calculation method (Field 1 of 
the PROBMETH keyword) is I, TOPO codes are: 

1 = ridgetop 

2 = sidehill 

3 = bottom. 

When the conditional probability calculation method is 2 or 3 (in Field 
1 of the PROBMETH keyword), TOPO codes are: 

1 = ridgetop or upper slope 

2 = midslope or lower slope. 

Field 1 : Topographic position code (default is 1). 

ASH DEPTH When the conditional probability calculation method (specified in Field 
1 of the PROBMETH keyword) is 2, the soil ash (loess) depth is needed 
to calculate the conditional probability that the subject stand will be in- 
volved in a regional outbreak. 

Field I : The ash depth in inches (default is 15.93). 

The following two examples illustrate the use of several of the preceding keywords in 
tailoring a simulation to specific situations. 


*Note that the available conditional probability of outbreak models were all developed using data describing 
stand conditions in only one regional outbreak— the 1971-74 outbreak. Such models predict conditional probabil- 
ity of outbreak only as a function of site and standi characteristics, even though the dependent variable ts also con- 
ditional on a number of unobserved factors, such as climate and weather. The fact that such climatic factors are 
very difficult to quantify and relate to specific outbreak histories (Mason and Luck 1978) does not make the prob- 
ability of stand outbreak models any less conditional on the specific climatic factors associated with (and perhaps 
peculiar to) the 197 1-74 outbreak. The potential for such models to overestimate the probability of a stand being 
Involved in future outbreaks is large, in our opinion. 
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Outbreak Initial 
Conditions 


Example i. 

You desire to simulate the following conditions. You hypothesize that the annual 
probability of a regional tussock moth outbreak occurring is 0.1, given that at least 20 
years has elapsed since the last tussock moth outbreak (which occurred in 1971). You 
want to use Heller’s model to stochastically determine whether or not there will ac- 
tually be an outbreak in the sample stand if there is a regional outbreak. You believe, 
however, that Heller’s model overestimates the conditional probability of stand out- 
break by a factor of 4. Furthermore, you are curious to see just what the “large 
amount of intermediate output detailing the operation of the Combined Model” 
looks like. Finally, you would like to use the default values of all other keywords. 

The following group of tussock moth keywords will accomplish this: 

DFTM 

RANSCHED 20. 0.1 1971 

RANSTART 

PROBMETH 1. 0 25 

DEBUG 

END 

Example 2, 

You would like to see how the results from the simulation described in example 1 are 
changed when you reseed the random number generator. You are also no longer 
curious to see the additional output DEBUG produces. Adding the RANNSEED 
keyword with any three odd numbers (and deleting DEBUG) will accomplish this: 

DFTM 


RANSCHED 

20. 

0.1 

1971. 

RANSTART 

PROBMETH 

1. 

0.25 


RANNSEED 

13. 

571. 

14327. 


END 


The DFTM Outbreak Model requires information describing tussock moth population 
levels and foliage biomass at the start of an outbreak (Phase I), on a 1000-in' midcrown 
sample branch basis. Recall that this sample branch represents a group or class of trees in 
the outbreak model. Larval density at the start of every outbreak can be assigned either 
randomly (RAN LARVA) or deterministically (DETLARVA) for each host species; the same 
method must be used for both host species. The BIOMASS keyword specifies how foliage 
biomass is to be determined. And if sample-based estimates of the mean and standard 
deviation of the foliage distribution are available, they can be incorporated via the DFBIO- 
MAS and GFBIOMAS keywords. 

RAN LARVA Randomly allocate first instar larvae to the tree classes at the start of 
every outbreak. Note that this is the default larval allocation method. 
The host species is specified in Field 1 ; larval density is drawn from 
a normal distribution with mean specified in Field 2 and standard 

deviation specified in Field 3. The mean larval density (Field 2) may 

vary from outbreak to outbreak if the between-outbreak standard 
deviation (Field 4) is positive. 

Field 1: Host species: 

1 = Douglas-fir, and 2 = grand fir. 

Field 2: 'The average number of first instar larvae per midcrown sample 

branch for the host species specified in Field 1; defaults are 9 and 11 

larvae per sample branch for Douglas-fir and grand fir, respectively. 
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Field 3: 

Field 4: 


DETLARVA 

Field 1: 

Field 2: 

Field 3: 

Field 4; 

BIOMASS 

Field 1: 


The within-outbreak standard deviation of first instar larvae; default 
is 2.0 larvae for both host species. 

The between-outbrcak standard deviation of first instar larvae; 
default is 0.0 for both host species. If this parameter is positive, the 
average larval density for the proper host species will be randomly 
chosen at the start of every outbreak by drawing from a normal 
distribution with mean specified in Field 2 and standard deviation 
specified in Field 4. Larvae will then be allocated to individual tree 
classes by randomly drawing from a normal distribution with this 
randomly chosen mean, and standard deviation specified in Field 3. 
This feature is best suited for use with the keywords that produce 
multiple outbreaks (primarily RANSCHED). 

Deterministically assign different levels or densities of first instar larvae 
to each third of the tree classes, after sorting the tree classes by average 
diameter (in descending order). 

Host species: 

1 = Douglas-fir, and 2 = grand fir. 

The number of larvae assigned to the largest third of the tree classes; 
defaults are 1 1 and 15 for Douglas-fir and grand fir, respectively. 

The number of larvae assigned to the middle third of the tree classes; 
defaults are 9 and 10 for Douglas-fir and grand fir, respectively. 

The number of larvae assigned to the smallest third of the tree 
classes; defaults are 7 and 5 for Douglas-fir and grand fir, 
respectively. 

Specify the method for calculating foliage biomass (in grams) and 
percentage new foliage on the 1000-in' midcrown sample branches; 
default is method 4. 

The calculation methods are: 

1 = The foliage biomass and percentage new foliage values will be 
randomly drawn from a normal distribution. The default mean and 
standard deviation of the foliage biomass distribution is 21 5g and 64g 
for Douglas-fir and 227g and 64g for grand fir, respectively (from 
Hatch and Mika 1978, p. 16). The default mean and standard 
deviation of the percentage new foliage distribution is 27 and 13 for 
Douglas-fir, and 35 and 7 for grand fir, respectively (from Hatch and 
Mika 1978, p. 21). These defaults can be replaced by using the 
DFBIOMAS and GFBIOMAS keywords. 

2 = The species-specific equations developed by Hatch and Mika 
(1978) are used deterministically to predict the percentage new 
foliage and foliage biomass from tree and site variables. Because the 
equations behave poorly, Monserud (i978a) recommended that this 
option not be used. 

3 = Species-specific equations developed by C. R. Hatch, Univer- 
sity of Idaho (see Monserud 1978a, p. 55) will be used to deter- 
ministically predict percentage new foliage and foliage biomass; the 
only independent variable is basal area percentile. 
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4 = Same as method 3, but with a random normal error (with 
mean ~ 0) added to each prediction. The standard deviation of this 
random error distribution equals the standard error of the regression 
fit described in method 3. For foliage biomass, the standard devia- 
tion is 57g for Douglas-fir and 58g for grand fir; for percentage new 
foliage, the standard deviation is 1 1 for Douglas-fir and 7 for grand 
fir (see Monserud 1978a, p. 55). 

DFBIOMAS 

Used to replace the default parameter values determining Douglas-fir 
sample branch foliage biomass, when method 1 is specified on the 
BIOMASS keyword; foliage biomass and percentage new foliage will 
be randomly drawn from a normal distribution with mean and standard 
deviation specified by Fields 1 through 4: 

Field 1: 

The mean of the foliage biomass distribution for Douglas-fir; default 
is 214g. 

Field 2: 

The standard deviation of the foliage biomass distribution for 
Douglas-fir; default is 64g. 

Field 3: 

The mean of the percentage new foliage distribution for Douglas-fir; 
default is 27. 

Field 4: 

The standard deviation of the percentage new foliage distribution 
for Douglas-fir; default is 13. 

GFBIOMAS 

Same as DFBIOMAS, but for grand fir; used in conjunction with 
BIOMASS method 1 . 

Field 1: 

The mean of the foliage biomass distribution for grand fir; default is 
227g. 

Field 2: 

The standard deviation of the foliage biomass distribution for grand 
fir; default is 64g. 

Field 3: 

The mean of the percentage new foliage distribution for grand fir; 
default is 35. 

Field 4: 

The standard deviation of the percentage new foliage distribution 
for grand fir; default is 07. 

Example 3. 

You would Uke to make a projection with only one DFTM outbreak. You also want this 
outbreak simulated in year 1971 (which happens to be when the first cycle of the projec- 
tion be^ns). You have estimated that the outbreak probably began in this stand with an 
average of nine first instar larvae per KXXhin' midcrown sample branch on Douglas-fir, 
and 12 larvae per s^ple branch on grand fir; your best estimate of the standard 
deviation is approximately four larvae per sample branch on either host species. You 
prefer to use biomass method 1 rather than the default method ( 4 ). The following key- 
words wUl accomplUh this (note that M ANSTART is supplied by default): 

DFTM 

MANSCHED 1 . 

RANLARVA 1 . 9 . 

RANURVA 2 . 12 * 4 ’ 

BIOMASS 1 . ’ ‘ 

END 



umpui f. . 1 I u * 

You want to modify the simulation in example 3 to include multiple outbreaks that are 
stochastically determined but occur approximately every 45 years, with a minimum 
waiting time of 36 years (note that these assumptions imply an annua! probability of 
regional outbreak of 0. 1). You have very little information on how outbreak severity 
varies in the long run, but you are sure that it is not constant; thus you assume that the 
between-outbreak standard deviation of larval density is approximately 5 larvae for both 
host species. You would also like to modify your foliage biomass assumptions as follows: 
Mean foliage biomass of a lOOO-in' midcrown sample branch is lOOg arid 2(X)g for 
Douglas-fir and grand fir, respectively, while mean percentage new foliage is 20 percent 
for Douglas-fir and 30 percent for grand fir; you have no information that warrants 
replacing the default foliage standard deviations. The following keywords will mimic 
these assumptions: 


DFTM 

RANSTART 


RANSCHED 

36. 

0.1 

1971. 

RAN LARVA 

1. 

9. 

4. 

RAN LARVA 

2. 

12. 

4. 

BIOMASS 

1. 


20. 

DFBIOMAS 

100. 


GFBIOMAS 

200. 


30. 


END 


Tree Class 
Compression and 
Redistribution Options 


Before the DFTM Outbreak Model is called, the list of up to 1,350 tree records carried by 
the Stand Prognosis Model is compressed into a maximum of 100 groups or classes of trees 
(see NUMCLASS and WEIGHT). The purpose of this compression is to save computer 
time, by combining trees that are similar (as far as the DFTM Outbreak Model is concerned) 
into the same tree class. Once these tree classes have been created, the rate at which in- 
sects are assumed to annually redistribute between tree classes can also be specified (see^ 

RED 1ST) It is anticipated that very few users will have need of the keywords discussed in 
this section (viz., NUMCLASS, WEIGHT, REDIST). The default values should be ade- 
quate for most applications. ... 

As previously mentioned, only two tree characteristics (for a given species) are important 
to the DFTM Outbreak Model: the percentage new foliage and foliage biomass of the mid- 
crown sample branch. Recall that these two foliage attributes are assigned to each tree using 
the procedures described with the BIOMASS keyword. The compression routine is simply a 
procedure for deciding which trees are most alike with respect to their foliage complements. 

Two keywords control this compression: NUMCLASS and WEIGHT, The relative 
importance of the two foliage characteristics is determined by WEIGHT. And the NUM- 
CLASS keyword determines the number of tree classes (Field 1 or 2, depending on the 

species) to be created, using the following procedure: , . c u 

t. The mean and standard deviation for each foliage characteristic are calculated. Each 
tree’s foliage characteristics are then “standardized” by subtracting the mean and dividing 

by the standard deviation of the appropriate foliage characteristic. 

2. A new attribute is then created for each tree; call this attribute A . This attribute is the 
weighted sum of the standardized foliage characteristics; the weights used to multiply each 
standardized foliage characteristic in this sum are specified on the WEIGHT keyword card. 
Attributed is then sorted into descending order, and used by the following two compres- 
sion algorithms to assign trees to tree classes: 

3. The first compression algorithm finds the largest gaps (or differences or distances) be- 
tween adjacent sorted A values. These gaps become the boundaries for compression. All 
trees between two adjacent gaps are classified or grouped into the same tree class. This 
algorithm is intended to find those trees (or groups of trees) that have foliage charactenstics 
that are very unusual, and should therefore not be combined with other trees into the same 
tree class. Generally speaking, this first compression algorithm does a good job of finding 
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DFTM Control and 
Stand Management 
Options 


such unusual (rees, but in the process creates a few tree classes that contain a large number 
of trees (which are relatively similar). Additional discussion of this algorithm is given by 
Monserud (1978a, p, 59-61; see rule 4B). The proportion of the total number of tree classes 
determined by this algorithm is specified by parameter 3 of the NUMCLASS keyword. 

4. The remaining available tree classes are determined by the second compression algo- 
rithm, which works as follows: The tree class containing the largest number of tree records 
is split into two classes, so that each class contains half of the records in the original class. 
Again the tree class containing the largest number of tree records is found, and then split 
evenly into two classes. This algorithm is repeated until the number of tree classes specified 
on the NUMCLASS keyword is created. The second compression alogrithm is intended to 
insure that one tree class does not contain an excessively large number of tree records, even 
though the foliage characteristics of those trees are relatively similar. 


NUMCLASS Specifies the number of tree classes to be created for each host species 
(Fields 1 and 2), and the proportion of tree classes to be created by the 
first tree compression algorithm (Field 3). All trees in a given tree class 
will be represented by one midcrown sample branch in the DFTM Out- 
break Model. 


Field 1 : Number of Douglas-fir tree classes; default is 20. 


Field 2: 


Field 3: 


WEIGHT 


Field 1: 
Field 2: 
REDIST 


Field 1: 


Number of grand fir tree classes; default is 20. (Note: the sum of 

Fields 1 and 2 must not exceed 1(X)) 

Proportion of tree classes determined by the first tree class compres- 
sion algorithm (see preceding discussion); default is 0.50. 

Specify the relative importance of the two foliage variables used by the 
algorithm for compressing the list of trees into the number of tree 
classes specified by the NUMCLASS keyword (see discussion preceding 
NUMCLASS keyword). The default values result in percentage new 
foliage and foliage biomass being equally important. 

The weight given to percentage new foliage; the default value is 1 .0. 

The weight given to foliage biomass; the default value is 1 .0. 

Specify the annual redistribution rate of insects between tree classes (see 
Colbert and Wong 1979, p. 54). In effect, the redistribution rate (Field 
1) operates by reducing the variation between the number of insects per 
tree class (weighted by the number of trees per tree class). The default 
redistribution rate of 0.25 will reduce the between tree class variation in 
insects by 25 percent for each year of the outbreak. A rate of 0.0 results 
in no redistribution, and a rate of 1 .0 results in completely uniform re- 
distribution, with the same number of insects in each tree class after the 
first year of the outbreak. 

The annual tussock moth redistribution rate; default is 0.25. 


The CH EM I CAL, N PV2, and NPV3 keywords are available for simulating the effect of 
applying either a chemical control or a virus, at various occasions during the outbreak. 

If the user desires to simulate a control measure that is not in the available list, then the 
TMPARMS keyword can be used to alter the appropriate mortality rates or growth param- 
eters in the DFTM Outbreak Model; in this case, Colbert and Wong ( 1979) must be con- 
sulted to calculate the appropriate parameter value. A SALVAGE option is also available. 
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CHEMICAL 


Chemical control will be applied in the phase of the outbreak specified 
in Field 1 to the instar specified in Field 2, and with the efficacy speci- 
fied in Field 3. Note that any number of CHEMICAL keywords can be 
used. 


Field 1: 


Phase (year) of the outbreak when chemical control will be applied; 
default is 3. 


Field 2; 
Field 3: 

NPV2 

NPV3 

SALVAGE 

Field 1: 

TMPARMS 


Instar that will be targeted for control; default is 4, 

Instar specific mortality rate resulting from the chemical control 
treatment; default is 0.95. 

Nuclear polyhedrosis virus will be applied in Phase II. 

Nuclear polyhedrosis virus will be applied in Phase III. 

At the end of a tussock moth outbreak, salvage all surviving host trees 
that have been defoliated more than the percentage tree defoliation 
specified in Field 1 . 

The minimum percentage tree defoliation for trees that will be 
salvaged; default is 50.0. 

This keyword has been provided for experienced users who have need 
to change additional parameters in the DFTM Outbreak Model. 

The parameters in question are most of those which make up the 
“PARAMETER” file as described in appendix A of Colbert and 
Wong (1979). The DFTM submodel described here contains an internal 
storage area which acts as a surrogate to their PARAMETER file; the 
TMPARMS keyword can be used to replace any parameter in this stor- 
age area. Colbert and Wong (1979) have defined the PARAMETERS 
and illustrated how new values are calculated. They have also prepared 
a table (see appendix A of their paper) which refers to the variables by 
data-card number, represented by the letter /, and value-on -the-car 
number, represented by the lettery. These same vdues, / andy. are used 
to reference the parameter values to be redefined m the Combined 
Model. Note that the parameter value will be reset to its default value i 
additional stands are processed in the same run. 

The card or record number (/, as described by Colbert and Wong 
1979 p 41-46), which contains the parameter to be replaced in the 
DFTM submodel storage area. The value of / must equal an integer 
from 2 to 12 or 19 to 25. 

The jth value on the ith card which corresponds to the ptuameler to 
be replaced in the DFTM submodel storage area. The value otj must 
equal an integer from 1 to 6. 

The value which is to replace the parameter corresponding to the Jth 

v^ue on dte fth card in the parameter fde. An error w, 11 occur tf thts 

or either of the preceding two fields are left blank. 

ToiUustratetheuseoftheTMPARMSkeywordta^lM 

TMPARMS keywords that or*NPV^s used the in^ar-specific daUy disease mor- 

found in Field 3 in table 1 . 


Field 1: 


Field 2: 


Field 3: 
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Table I.~Parameter values on the TMPARMS keyword(s) that will mimic the virus control 
keywords; note that four TMPARMS keywords are needed to mimic each of the 
virus control keywords 


Keyword 
to be 

Parameters of the TMPARMS keyword 

mimicked 

Field 1 

Field 2 

Field 3 

NPV2 

3 

3 

0.036 


3 

4 

.039 


3 

5 

.042 


3 

6 

.072 

NPV3 

4 

3 

.036 


4 

4 

.039 


4 

5 

.042 


4 

6 

.072 


Example S, 

You would like to rerun the simulation in example 3, but with the following addi- 
tions: simulate a chemical control with 90 percent efficacy applied in Phase III of the 


outbreak to the second instar, and salvage 
percent: 


DFTM 

MANSCHED 1 . 

RANLARVA 1 . 

RANLARVA 2. 

BiOMASS 1 . 

CHEMICAL 3. 

SALVAGE 75. 

END 


II trees that were defoliated more than 75 

9. 4. 

12. 4. 

2. 0.90 


INFORMATION PRODUCED 


The Combined Model displays a variety of output tables which summarize the opera- 
Uon of the DFTM Outbreak Model during the course of the simulation. The DFTM 
Options and Input Table summarizes and describes the keywords that are in effect dur- 
mg t e simulation. The DFTM Outbreak Summary Table lists the information germane 
to the scheduling and timing of outbreaks. The DFTM Defoliation Statistics Table— 

w 1 C IS pro uced for each outbreak— displays the effect of tussock moth defoliation on 
each tree class. 


I hew tussock moth outputs supplement the normal output tables produced by the 
Stand Prognosis Model (see Stage 1973, and Wykoff et al. 1982 for examples). An ex- 
^ple of each t^e of output wiU be provided in the following discussion. A listing of 
e runs ream that produced the simulation output illustrated in this section can be 
tound in appendix A. 


DFTM Options and 
Input Table 


in ^ sy"^^ary of the important tussock moth keywords and parameter values used 

« contained in the DFTM Options and Input Table (fig. 3). 
and ^ wc listed m the left-hand column. A short description of the keyword 

the numenc parameter values used in the simulation then foUow to the right, com- 
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DOUGLAS-FIR TUSSOCK MOTH IN DOUGL^S-FIR AND GRAND FIR: DFTM VERSION 3.1; PROGNOSIS (INLAND EMPIRE) 
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Hgure 3. — Sample output from the Combined Model: DFTM Options and Input Table. 



rri;ing the main body of (he (able. Note that the DFTM Options and Input I’ahic ix 
de^ened to display only the options in effect for the simulation; it doe.s not include 
sliotjlie available keywords. Also note that the DFTM keywords that arc explicitly 
specified by the user are also listed (in the order specified) in the “Options Selected by 
Input” table in the Stand Prognosis Model (see Wykoff et al. 1982 for an example). 




The Df-TM Oulbreak Summary Table (Hg. 4) displays the timing of liissoek iiiotli 

muba-aks in summary form. The Combined Model contains logic tlmt controls tlic 

f J""® periods) will contain a tussock motli 

< ‘"'P '‘"P® ‘P "lis decisio .1 proces,s, Tlic nr.sl stun 

re>'icinai outhreT"®'?! and the second detennine.s which of tlie 

rCfc-onaJ outbreaks will include the subject stand. 



- DFTM OUTBREAK 

SUMMARY TABl E 


; CYCLE.' 

‘‘■•''••f’'' YEARS 

YEAR OF 
REGIONAL DFTM 
OUTBREAK 

CONDITIONAL 

PROBAB 1 1. 1 TY 

OF STAND OUTBREAK 

WAS flfCRf' 

<3i/fbri:ai< 

STAND YBIf>- 

..' ^571 - 1976 

■■ H'/C - 1981 

19.J1 - 1991 
- 1<;91 - P001 

> 2i)!,n - 2006 

ri d' 0 f !6 - 20 n 

' i:0n - 2021 

1971 

2006 

0.81*3 

0 . 1*96 

0.5112 

0. 661* 

0 . 71*0 

0.763 

0.1*18 

YES 

NO 

NO 

NO 

NO 

YES 

NO 


Afi 

it{ 

I :'3? 


t'i!rt.re4,- 

TaM?. 


“odd: DFTM Outbreak Summary 


the outcome of the first step of the outbreak Summary Table (fig. 4) list 
MiKlel cycles is contained in the first two enh, ^ Prognosis 

'■^?‘onaI outbreak years. When the RANsSdPn column contains a list of 

wniple, a senes of regional outbreak dates 

da es are assigned to existing proje fon cvlT generated. These outbreak 

iS searsof ,he regroup ouIbre~rwl a '^iri.it. 

The remainder Of the table summari, ® inserted. 

■nduded'-'’J'h''"’" ^ °f conditfond hT°'^“’u ricci.sion logic, 

cveie ev.n lu ' ^'gtonal outbreak. Note that a cnnrfc^^' 1^'** subject stand will be 

■ ' though there may be no reeinnai P^ditional probability is printed for every 

ZtjUr P^babSe us“ed "t I" 

thTS;' 7 ■" ‘.he example presented he e) Thets? ^^NSTART option is 

“ —1 1. ...Ut 



Mode?,o sim“|a"f 

•> »™,„g mesage will be pri„,™' “ lung tussock moth outbreak. 


On 



DFTM Defoliation 
Statistics Table 


Figure 5 illustrates the output summarizing each tussock moth outbreak . The first sever; 
lines contain a message indicating which years of the subject stand’s development include 
the outbreak that is detailed in the rest of the table. A message is also printed indicating 
when a salvage is scheduled, 

The main body of figure 5 consists of the ' ‘Summary of Tree Class Characteristics.” A 
row of summary statistics is printed for each tree class, which is a collection of trees with 
similar foliage attributes— the only tree characteristics important to the DFTIVI Outbreak 
Model, These statistics are divided Into two groups: “Before Outbreak” (columns l-IO) 
and “After Outbreak” (columns 1 1-19). The last three rows contain weighted averages of 
the same characteristics for each host species (Douglas-fir and grand fir) and for all host 
species combined. Except for the columns labeled “Records per tree class” and “Trees pei 
acre” (columns 3 and 4), all summary statistics printed in this figure are averages weighted 
by the number of trees per acre (column 4) represented by each tree in a given tree class. 

Percentage branch defoliation (column 1 1) is the only variable in figure 5 directly pre- 
dicted by the DFTM Outbreak Model. Percentage branch defoliation is converted to per- 
centage tree defoliation (column 12) using the function illustrated in figure 2. The insensith 
ity of percentage tree defoliation to any amount of branch defoliation below 55 percent car 
be seen from examining columns 1 1 and 12 in the “Summary of Tree Class Characteristics’' 
table illustrated in figure 5; the extreme sensitivity of the relationship graphed in figure 2 is 
also appai ent for values of branch defoliation between 60 and 90 percent. The predictions 
of mortality, top-kill, and diameter and height growth loss are all functions of percentage 
tree defoliation. 

Recall that tree mortality (or survival) is modeled as a continuous ratlier than a discrete 
event In the Stand Prognosis Model (Stage 1973; Monserud 1978a). Mortality operates by 
reducing the number of trees per acre represented by a given tree record. Thus the number 
of trees per acre in a given tree class after the outbreak equals the product of the trees per 
acre before the outbreak (column 4) and 1.0 minus the 5-year mortality rate (column 13). 
Of course it would not be correct to attribute the reduction in trees per acre solely to the 
tussock moth, for the normal mortality rate in the absence of a tussock moth outbreak is 
certainly greater than zero. 

Top-kill is not modeled in the same manner as mortality however, for top-kill is treated ai 
a discrete event in the combined model. Monserud (1978a, p. 67-68) detailed the procedure 
for stochastically determining whether or not a given tree will have top-kill, and how much. 
The average amount of top-kill for each tree class is summarized in columns 18 and 19 of 
figure 5. As a result of top-kill, the total height, live crown ratio, and volume of a tree are 
reduced accordingly; the procedure for calculating the volume of a top-killed tree is dis- 
cussed by Monserud (1980, 1981). 

The ef fect of tussock moth defoliation on diameter and height growth can be seen in 
columns 14-17 of figure 5; the sum of net growth and growth loss equals the growth in the 
absence of tussock moth defoliation. Note that the change in height columns (16 and 17) 
include top-kill losses. Because top-kill losses can potentially exceed heiglit growth, net 
change in height can be negative (see column 16 in figure 5 for tree classes 9 and 22). The 
average amount of top -kill as a percentage of live crown length is displayed in column 19 of 
figure 5. 
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Kher Output 


As discussed earlier, the defoliation estimates made by the DFTM submodel are auto- 
matically transmitted to the Prognosis Model and translated into tree damage. Thus, any 
tree mortality, growth reductions, and volume losses attributable to tussock moth are re- 
flected in the normal output produced by the Stand Prognosis Model. 

Figures 6 and 7 are examples of the two major displays produced by the Prognosis Mo< 
these output tables are discussed by Stage (1973, p. 3-5) and Wykoff et al. (1982). An exj 
ination of figure 6 reveals that mortality losses were estimated to be 400 ft’/acre/yr durini 
the simulated outbreak that began in 1971, During this outbreak, 52 percent of the volim 
lost was grand fir and 45 percent was Douglas-fir (see the species composition summary o 
the right half of figure 6). The outbreak was followed by a salvage that removed 989 ftVa 
in the year 1976 (all trees with tree defoliation exceeding 90 percent were .salvaged). In figi 
7 it can be seen that stand basal area was reduced from 173 ftVacre before the 1971 out- 
break to 99 ft Vacre after the outbreak; the salvage cut reduced it further to 57 ftVacre. I 
evident from both figures 6 and 7 that the combined model predicted that the stand con- 
tained 618 trees per acre before the 1971 outbreak, 382 trees per acre after the outbreak, a 
210 trees per acre after the salvage at the end of the 1971 outbreak. 

Finally, optional tussock moth output is printed if the user specifies the proper key- ; 
word(s); see the descriptions of the DEBUG, DEBUT REE, and PUNCH keywords for 
details. 
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BEHAVIOR OF THE COMBINED MODEL 


The purpose of this section is to illustrate the behavior of the Combined Model when 
the major input options are varied. This section does not purport to be a sensitivity 
analysis, however. In all figures in this section, simulated stand volume development is 
plotted against time over the course of a 50- to 100-year projection, using the same 
stand. The sample stand simulated is B. E. Wickman’s Y-Ridge plots 1, 2, and 3, as 
measured at the beginning of the 1971-74 Oregon Blue Mountains outbreak. This stand 
is two-storied with 70 percent of the trees less than 6.8 inches d.b.h. in 1971 (before the 
outbreak). The two host species for tussock moth comprise most of the stand: grand fir 
and Douglas- fir account for 72 and 21 percent of the trees per acre and 46 and 29 per- 
cent of the total volume in 1971, respectively, 

The “no outbreak” curve will be the same on all figures, and indicates what the expected 
volume development over time would be for this stand if an outbreak had not occurred. 
Since the version of the Stand Prognosis Model used in this analysis did not include regener- 
ation establishment, all simulation results reported here are based on projecting only trees 
and seedlings already established at the start of the outbreak. Thus, the simulations that 
produced heavy mortality or large salvage removals quite likely underestimate stand volume 
near the end of the 50-year projections. 

Keep in mind that these examples are projections into the future, which obviously in- 
vo ves considerable uncertainty; our best guess should never be confused with certainty. 

so keep m mind that this one sample stand (or any one stand) is unlikely to be typical of 
an ^ea as large and diverse as the Blue Mountains of Oregon and Washington . 
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After 45 years, projected volume development following the lowest two levels of out- 
break (5 and 8 larvae per tree class) illustrated in figure 8 nearly caught up with the no- 
outbreak volume level. However, even 50 years after the two most severe simulated out- 
breaks began (14 and 20 larvae), volume development was still 700 to 1300 ftVacre behind 
the no-outbreak level. Obviously, insect density at the start of the outbreak can influence 
the simulated stand’s development far into the projection. 

This figure illustrates an additional important point. Consider one of the most severe out- 
breaks graphed, the 14 larvae per tree class curve. Five years after the outbreak began, 
almost 2{X)0 ftVacre were estimated to be lost. But 50 years after the outbreak began, only 
700 ftVacre were lost (and probably less than that since any regeneration subsequent to the 
outbreak was ignored). If this particular stand were not scheduled for harvest for, say, 50 
years, then it is misleading to tell the manager that 2000 ftVacre have been lost, for the 
stand will probably “find” much more than half of these lost cubic feet if left alone for 45 
more years. The point of this admittedly oversimplified argument is that “loss” is not cons- 
tant over time. Perhaps more realistic measures of loss would be either the expected number 
of additional years required to reach the volume anticipated at the end of the rotation if 
there had been no tussock moth outbreak, or the expected volume loss at the end of the 
scheduled rotation. 



Figure 8.— Comparison of simulated volume development resulting from outbreaks 
beginning with an average of 5, 8, 14, 20, and 100 first instar larvae per tree class. 
Larvae were randomly allocated (using RAN LARVA) with a within-outbreak stand- 
ard deviation of 2 larvae. 
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Figure 10.— Simulated volume development resulting from outbreaks beginning with 
an average of 5 first instar larvae per tree class. Larvae were randomly allocated with a 
standard deviation of 0, 2, 4, 8, and 16 larvae. 
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Figure 1 L— Simulated volume development following outbreaks resulting from deter- 
ministically allocating three larval densities (Douglas-fir: 7, 9, 11; grand fir: 5 , 10, 15) 
to each third of the average diameter distribution. Each outbreak resulted from the 
same total number of first instar larvae. Note that each tree class represented an un- 
equal number of trees per acre , 




All possible combinations of allocating these larval densities to the three diameter classes 
were simulated to produce figure 1 1 . The most severe reduction in volume was obtained by 
assigning the high larval density to the middle third rather than the top third of the diameter 
distribution. This apparent anomaly occurred because the middle third of the diameter ^ 
distribution represents more trees per acre (and more volume) than the upper third. In light 
of this, it is not surprising that the least severe reductions in volume resulted from assigning 
the low larval density to the middle third of the diameter distribution. Recall that the 
number of trees per acre each tree class represents is determined both by the original stand 
inventory design and the method for allocating tree records to tree classes (Field 3 on the 
NUMCLASS kej^vord record); the methods for allocating lai*vae to the tree classes are in- 
dependent of the number of trees per tree class. 



Figure 12.— Simulated volume development following outbreaks resulting from deter- 
ministically allocating three larval densities (Douglas-fir; 7, 9, 11; grand fir: 5, 10, 15) 
to each third of the average diameter distribution; each tree class represents approx- 
imately an equal number of trees per acre. Each outbreak resulted from the same total 
number of first instar larvae. 

In the simulations illustrated in figure 12, larvae were allocated deterministically to tree 
classes exactly as they were in figure 11. For this case, however, biomass was allocated de- 
terministically (using BIOMASS method 3 rather than the default method 4), and the 
successive-halving rule was used exclusively to create tree classes (Field 3 was 0.0 on the 
NUMCLASS keyword). These two changes resulted in an approximately equal number of 
tree records in each tree class, and thus the number of trees per acre in each third of the 
diameter distribution was roughly equal. The result of varying the allocation of a fixed 
number of larvae by tree size-class can be quite dramatic, as can be seen in figure 12. As ex- 
pected, the niost severe outbreaks (judging from volume loss) occur when the high level of 
larvae are assi^ed to the largest third of the diameter class, and the least severe outbreaks 
occur when this high level of insects is distributed to the lowest third of the diameter class. 

This latter case is quite interesting, for the standing volume 45 years after the lightest out- 

® outbreak’’ simulation, even though 

t /acre were killed by the tussock moth in this outbreak. Concentrating the larvae in 
the smaller diameter trees is apparently silviculturally similar (in this case) to a light thinning 
m a stand that is overstocked in the small diameter classes. Although tussock moth out- 
reaks mvanably result in a short-term volume loss, these results indicate that some out- 
breaks have the potential for mimicking wise management, and being beneficial in the long 



run. It is quite apparent from figure 12 that stand volume development can vary greatly, 
depending on the size of tree defoliated during the outbreak. Using overall stand averages 
for insect densities thus results in less precise estimates of volume development if the pest is 
differentiating between tree size classes, or redistributing more to trees of one size class than 
another. 

The effect of varying the method of allocating foliage biomass to the tree classes is ex- 
amined next (fig. 13). In this and subsequent runs, the random larvae option was used to 
allocate an average of 14 first instar larvae per tree class, with a standard deviation of 2 lar- 
vae. A fairly large amount of variability results from varying the foliage biomass method. 
This should not be surprising, for foliage biomass is the only tree characteristic important to 
the DFl'M Outbreak Model. Note, however, that all the outbreak volume-over-time curves 
in figure 13 converge, except that which was produced by using BIOMASS method 2 (curve 
B); this method uses the deterministic biomass equations of Hatch and Mika (1978) that 
Monserud (1978a) concluded are poorly behaved. 
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Figure 13,— Comparison of simulated volume development resulting from outbreaks 
that began with each of the four methods of allocating foliage to the lOOO-m mid- 
crown sample branch that represents a tree class. In, this and all subsequent figures, 
the random allocation method was used to distribute an average of 14 first Instar lar- 
vae per tree class, with a within-ouibreak standard deviation of 2 larvae, except for 
figures 22 and 24. In addition, the default foliage biomass option (method 4) was used 
for all other figures. 


Figure 14 illustrates the effectiveness of simulating various control measures for the same 
outbreak conditions. The ‘ ‘no control ‘ curve (G) in figure 14 is the same as curve E in 
figure 8 and curve D in figure 13: foliage was determined by the default biomass option 
(method 4), and first instar larvae were assigned to tree classes randomly (average = 14, 
standard deviation = 2). Control treatments (especially virus) that are applied in Phase II of 
this simulated outbreak almost completely control the tussock moth (as far as volume loss 
is concerned), unless the treatment is of low efficacy (e.g. , 80 percent) and applied late m 



Phase II (curve C). Delaying control treatments until Phase III results in much greater varia- 
tion in effect iveness. An early application of chemical control with a high efficacy (e.g., 95 
percent, curve B) is quite effective in reducing volume loss, while delaying application until 
late in the outbreak results in almost no control even through efficacy is high (curve F). In- 
termediate levels of control are obtained by applying virus (curve E) or chemical controls of 
lower efficacy (curve D) early in Phase III. All simulated control measures had one obvious 
effect in common: the later the application and/or the lower the efficacy, then the greater 
the volume loss at the end of tlie outbreak . 

The effect of varying salvage intensities is next shown (fig. 15). Note that the “no sal- 
vage’' outbreak curve (B) is the same as the “no control’’ curve (G) in figure 14. The major 
effect of salvage appears to be a reduction in the rate of stand development; the rate of 
reduction increases as salvage intensity Increases. Keep in mind that the volume available at 
the end of these 50-year projections is most likely underestimated for the heavy salvage 
treatments, because the current version of the Stand Prognosis Model does not yet simulate 
the development of seedlings that would become established subsequent to the salvage. It 
is apparent, however, that the time required for a simulated stand to reach preoutbreak 
volume levels can be lengthened considerably by increasing the salvage intensity. 

In figure 16 the annual insect redistribution rate (Field 1 on the REDIST keyword record) 
is varied. As in previous figures, first instar larvae at the start of the outbreak were allocated 
to tree classes randomly with a mean of 14 and a standard deviation of 2 larvae per tree 
class. As expected, the severity of the outbreak increases as the insect redistribution rate is 
increased from 0.0 (no annual insect redistribution) to 1 .0 (completely uniform redistribu- 
tion of insects among tree classes). Although not illustrated, varying the redistribution rate 
with a higher standard deviation (namely 8) for allocating larvae to tree classes produced 
essentially the same results as exhibited in figure 16. 



Figure 14.— Simulalcd volume development following application of various DFTM 
control measures. Virus and chemical controls applied in Phase II or III were simu- 
lated. In addition, two efficacy levels (80 percent and 95 percent) and early (instar 2) 
and late (instar 5) application of chemical control were considered. 




Figure 15.“Simulated volume development following various salvage intensilies. Sur- 
viving trees with defoliation exceeding 95 percent, 90 percent, 80 percent, and 50 per- 
cent were salvaged, Note that an additional 1835 ftVacre in mortality is also available 
for salvage in 1976. 



Figure 16.— -Simulated volume development resulting from varying the annual insect 
redistribution rate during the outbreak. A rate of 0.0 results in no redistribution, and 
a rate of 1.0 results In perfectly uniform redistribution of insects among tree classes; 
the default value is 0.25. 



In figures 17 though 20, the factors that affect the compression of the list of trees into 
tree classes are varied. The WEIGHT keyword is used to specify the relative importance of 
the two foliage biomass variables (i.e. , percentage new foliage and total foliage biomass 
of the sample branch). Curve C in figure 17 was produced by the default values for the 
WEIGHT keyword parameters: percentage new foliage and total foliage biomass both had a 
weight of 1.0 (i.e., percentage new foliage was considered just as important as total foliage 
biomass). Most combinations of WEIGHT keyword parameters produced simulated results 
that are bracketed by the two extremes: giving percentage new foliage no weight (curve B), 
and giving total foliage biomass no weight (curve F) in the tree class compression routine. 
An examination of figure 17 reveals that varying the parameters of the WEIGHT keyword 
resulted in minor changes in predicted volume for this particular stand. This is primarily a 
result of the method of allocating the foliage biomass to the trees (method 4 in this case). As 
long as there is a strong correlation between percentage foliage biomass and total foliage 
biomass, then changing the parameters on the WEIGHT keyword should have a minor ef- 
fect on the simulated outbreak. 



Figure 17.~SimuIated volume development resulting from varying the parameters on 
the WEIGHT keyword that detennine the importance of percentage new foliage rela- 
tive to total foliage biomass in the tree class compression algorithm. 

The number of tree classes (the first two parameters of the NUMCLASS keyword) used 
to represent the list of trees in the DFTM Outbreak Model is next varied (fig. 18). Observe 
that the default values of 20 tree classes per host species (curve C) result in a volume- over- 
time curve that is quite close to the curve produced by using 50 tree classes per species (curve 
D). Based on numerous projections, 20 tree classes per species does appear to provide a 
good approximation to the actual tree list. Note that using 5 tree classes per species (curve 
B) resulted in a volume-over-time curve almost coincident with the default 20 tree class 
per species curve (C). It would be incorrect, however, to infer that using 5 tree classes per 
species is as accurate as using 20 (although it is true in the example summarized by figure 
18). Generally, using less than 10 tree classes per species results in quite erratic behavior, and 
is not recommended. 




Figure 18.— Simulated volume development resulting from varying the number of tree 
classes into which the list of trees is compressed before calling the DFTM Outbreak 
Model. Half of the tree classes were created using the (first) maximum difference 
algorithm, and the remainder were created by the (second) successive halving com- 
pression algorithm. 


The third parameter (temporarily call it PROP) on the NUMCLASS keyword record 
specifies the proportion of tree classes to be determined by the first tree class compression 
algorithm (the majximum difference algorithm); the remainder are created using the second 
(successive halving) algorithm. To see the effect of varying PROP as the number of tree 
classes per species is also varied, compare figure 19 (using PROP = 0: the successive halving 
algorithm) and figure 20 (using PROP = 1: the maximum difference algorithm) with figure 
18 (PROP = 0.5, the default). It is apparent that exclusive use of the maximum difference 
algorithm (fig. 20) results in large variation in the consequent simulated outbreak as the 
number of tree classes is varied; the volume-over-time curve does not stabilize until the 
number of tree classes gets close to 50 per species. In contrast, the second compression 
algorithm (used exclusively in figure 19) produces stable results with much fewer tree classes 
(and less computing cost). 

Results obtained from simulated outbreaks that utilize the pseudorandom number gener- 
ator (Marsaglia and Bray 1968) are obviously conditional upon the sequence of random 
numbers used in the various calculations. The amount of variability in a given simulation 
that is due to the random number sequence is not at all obvious, however. By reseeding a 
given simulation several times (with the RAKNSEED keyword) and holding aU other input 

conditions constant, this “random” variability can be isolated and assessed. 
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Figure 20.— Simulated volume development resulting from varying the number of tree 
classes to be created by exclusively using the maximum difference algorithm (param- 
eter 3 on the NUMCLASS keyword equals 1 .0). 


An indication of the variability associated with the random larval allocation option is il- 
lustrated in figures 21 and 22. In figure 21, an average of 14 first instar larvae were randomly 
allocated to each tree class, with a standard deviation of 2 larvae; in figure 22 the standard 
deviation was increased to 8 larvae. Only the seeds of the random number generator were 
varied for the groups of simulations displayed in figures 21 and 22. As expected, the vari- 
ability in simulated volume increases with the standard deviation of the distribution that the 
random number generator is being used to produce. Keep in mind that the default biomass 
option (4) contains an additional source of random variation in figures 21 and 22, as does 
the algorithm for assigning top-kill damage to individual trees. 

When working with a simulation model containing varying degrees of random variabil- 
ity, there is always a risk that one simulation from a given set of initial conditions may be 
atypical. A good way to deal with this random variability is to generate replicate simulations 
(varying only the random number seeds) until the mean of the variable of interest (total 
volume over time in the examples in this section) is determined with acceptable accuracy; 
the user must of course consider the trade-off between the increase in both cost and preci- 
sion to be obtained from additional replicate simulations in determining what level of ac- 
curacy is acceptable. An even better way to deal with the variability in the simulated system 
is to consider both the estimated mean and variance of the variable of interest in decision- 
making. For example, output from replicate runs of the simulation model can be used to 
estimate the probability that a given critical value (say, 4,000 ftVacre lost) will be exceeded 
in deciding whether or not to apply a control measure. 



Figure 21.— Simulated volume development resulting from varying the sequence of 
random numbers used by the tussock moth related routines in the combined model. 
First instar larvae were allocated randomly to tree classes with an average of 14 and a 
standard deviation of 2 larvae. 
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Figure 22. — Simulated volume development resulting from varying the sequence of 
random numbers used by the tussock moth related routines in the combined model. 
First instar larvae were allocated randomly to tree classes with an average of 14 and a 
standard deviation of 8 larvae. 


The portion of the Combined Model where unpredictable variation is most evident is in 
the algorithms that stochastically control the timing and occurrence of outbreaks (see the 
RANSCHED, RANSTART, and PROBMETH keywords). To illustrate (fig, 23), consider 
the following simulation. Regional outbreaks were scheduled with a probability of occur- 
rence of 0. 1 , given that a 30-year waiting time had elapsed since the last regional outbreak. 
Whether or not the subject stand would be included in a regional outbreak was also deter- 
mined stochastically, based on the stand’s susceptibility to tussock moth. First instar larvae 
were allocated randomly to the tree classes (average = 14, standard deviation = 2); note 
that long-term larval density was held constant (between-outbreak standard deviation 
- 0,0 larvae). Default values were used for all other keywords and parameters. Five simula- 
tions were run; the first (fig. 23, curve B) used the default seeds for the random number 
generator, and the next four (fig. 23, curves C-F) used different random number generator 
seeds. 

It is clear from figure 23 that the variability associated with this stochastic outbreak op- 
tion is immense. Furthermore, it is quite unlikely that any given simulation using these same 
initial conditions could be called t>T)ical, If one is interested in estimating expected volume 
development over time resulting from these initial conditions, then it is almost essential to 
calculate the “average outbreak” volume curve (G). 

^ Figure 23 can also be examined in a different light. Recall that the large variability ex- 
hibited by curves B-F results from changing only the sequence of random numbers. It 
follows that any manipulation of the stand — no matter how minor — that in effect alters this 
sequence of random numbers may produce results comparable to reseeding the random 
number generator, Thus the projected yields for two stands that are very similar — but not 
identical — could be quite different if severe but infrequent stochastic outbreaks (such as 
in figure 23) are being simulated. The importance of replicating a given simulation with dif- 
ferent random number seeds cannot be overemphasized, especially when using the sto- 
chastic outbreak feature. 




Figure 23 —Five replicates ofsiimilniccl voKnnc devcInpmeiH rosuliiiig from viiryiiiB 

(he sequence of rniKlom luiinbers iisccl l>y iltc ins^uck nioiludiucd routines ill the 
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Figure 24 illusf rales six sucif average oui break volume curves, willi bolh the frequency 
and severity of the outbreaks being varietl. As in Urc previous figure, 
outbreaks was st ochastically determined by insing the RANSCH ED and RANSTART 
keywords, and larvae were randomly allocated lo tree classes. Three levels of seventy were 
simulated: an average of 4, 8. and 14 Urst insiarlurvacper tree class were allocated with a 
within-outbreak standard deviation of 2 larvae {roughly corresponding to light, modtmic, 
and heavy severity, respectively). As In die previous Hgiirc, long term larval density w^ held 
constant (Field 4 on the RANLARVA key word record wasO.O). For each 
of severity, both frequent and infrcquciu outbreaks were scheduled: the RANSUHtu 
parameters for the frequent outbreaks were 7 years minimum time between outbreaks with 
a subsequent annual probability of regional outbreak of 0.3; the correpondmg parameters 
for the infrequent outbreak schedule were 30 years and 0. 1 , respectively. 
value of the time between outbrcak.s was thus 9-2/3 years for the frequent out brea - 

lations and 39 years for the infrequent outbreaks. For each of the sixcoinbinations o o 
break severity and frequency, five replicates were simulated (by varying the seeds of the 
random number generator). Note that each curve graphed in figure 24 is le averag 
these five replicates. The default values for all other keywords were used. 

For the range of conditions examined in this example, average volume os ( 
sock moth) over time increased as either the severity (he., the number of larvae a 
of an outbreak) or the frequency of the simulated outbreak was increase , is is a 
predictable result. What is not very predictable, however, is the . .j. 

quency versus severity of outbreaks on long-term volume yields. It is quite - „ . 
generalizations in this regard, for figure 24 contains both examples of _ 

breaks being a more important factor In explaining volume loss than the number of larvae 
the start of an outbreak (compare cuves E and F) and vice versa (compare curves 




Figure 24.— Average simulated volume development resulting from varying the fre- 
quency of outbreaks for three different larval densities. Each outbreak curve (B-G) is 
the average of five simulations resulting from different random number sequences. 

The between-out break larval standard deviation was zero in these simulations. 

Furthermore, there are probably many combinations of these factors that will produce the 
same amount of volume at a given age or time. For example, curves C and D are quite close 
together, even though outbreaks in curve D began with twice as many larvae as those in 
curve C, and the average interval between outbreaks is 30 years longer in curve D than in 
curve C. Such results have implications for hazard or risk rating systems. An accurate 
hazard rating scheme should obviously consider the long-term probability of outbreak as 
well as the likely severity of individual outbreaks. It is the joint effect of frequency and 
severity of outbreaks that determines the expected loss over time due to tussock moth. 

The amount of variability associated with each of the average outbreak curves graphed in 
figure 24 decreased if either the severity of the simulated outbreaks was decreased or the fre- 
quency of the outbreaks was increased. Thus the average outbreak curve in figure 24 with 
the most variability is F; note that this curve is the same as curve G in figure 23. 

An important point can be made regarding the simulations summarized in figure 24: a 
few of the curves are likely unrealistic. For example, it is highly unlikely that the combina- 
tion of severe and frequent outbreaks averaged to produce curve G (and probably curve E) 
would ever occur in the same stand, although outbreaks have occurred as frequently in the 
Palouse Range and more severely in the Blue Mountains. Even though the same stand may 
be involved in repeated outbreaks, it is unlikely that such outbreaks would be equally severe 
and repeatedly start with the smne larval density. The random larval allocation method con- 
tains an option that allows average larval density to vary from outbreak to outbreak. This 
option was used to produce figure 25: the between-outbreak standard deviation was set at 
6,0 larvae (Field 4 on the RAN LARVA keyword record) and all the simulations that were 
used to produce figure 24 were rerun. Thus figure 25 illustrates six average outbreak curves 
(again based on five replicates per curve) with both within- and between- outbreak severity 

(i.e., larval density) as well as outbreak frequeneyvarying, 

The maj or difference between figures 24 and 25 is that outbreak frequency is far more im- 
portant when the between-outbreak larval standard deviation is moderately large (namely, 

6 larvae in fig. 25), All three of the f Trequept' ’ outbreak curves in figure 25 (namely C, E, 
and G, with an expected interval between outbreaks of 9,7 years) were below the “infre- 
quent outbreak curves (B, D, and F, with 39 years the expected interval between out- 




Figure 25.— Average simulated volume development resulting from varying both the 
frequency and the severity of the outbreaks. All simulations that produced figure 24 
were rerun after raising the be tween-outbreak standard deviation to 6.0 larvae (Field 4 
on the RAN LARVA keyword record). Each outbreak curve (B-G) is the average of 
five simulations resulting from different random number sequences. 

breaks). It may seem surprising that even curve C (frequent outbreaks with mean larvae 
= 4) is below curve F (infrequent outbreaks with mean larvae = 14); the reason is straight- 
forward, however. With frequent outbreaks occurring four times as often as infrequent 
outbreaks, the chances of getting a severe outbreak (e.g., mean larval densit y = 14) are 
greater than 1 in 4 when the between-outbreak standard deviation equals 6 larvae, even 
though long-term mean density is only 4 larvae . 

A final point suggested by both figures 24 and 25 concerns the use of expected yields (i.e., 
stand volume over time) in forest management planning. In an area where tussock moth has 
historically been a factor affecting stand development, it is probably overly optimistic to 
anticipate volume yields indicated by the ‘ 'no outbreak’ ’ curve. Although it may be diffi- 
cult to state with confidence which “outbreak” curve is most likely, almost all possibilities 
will predict less future volume than the “no outbreak” curve. The Combined Stand Prog- 
nosis/DFTM Outbreak Model has potential to reduce this bias associated with projecting 
future volume yields in stands susceptible to tussock moth. 

The features of the Combined Model that allow both the frequency and severity of out- 
breaks to vary stochastically in multiple-outbreak simulations were added with hopes of 
making long-term projections more realistic. Unfortunately, there is little information in- 
dicating how the severity (i.e., number of larvae) of an outbreak is distributed over time, 
just as there is little— if any— information available describing the probability of stand out- 
break for many different outbreak periods. Because of this dearth of knowledge, it was our 
objective to give the user considerable flexibility in stating assumptions regarding long-term 
interactions between the tussock moth and the stand being managed. The ultimate solution, 
of course, would be the development of a tussock moth population model that would b e 
truly dynamic; this would eliminate the need for the Combined Model to predict the prob- 
ability of outbreak, for the construct “outbreak” would then become an unnecessary 
artifact. 




The numerous simulations graphed in this section were intended to serve merely as 
examples of model behavior rather than as definitive statements that will hold in all cases or 
for all stands, although some generalization from a number of the figures would be war- 
ranted. Our hope is that the potential user will have a greater appreciation for the scope 
and versatility of the Combined Model after viewing the numerous scenarios that can be 
simulated. 
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